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ABSTRACT 

We analyzed the full Stokes spectra using simultaneous measurements of the photospheric (Fe I 630.15 
and 630.25 nm) and chromospheric (Mg I hi 517.27 nm) lines. The data were obtained with the HAO/NSO 
Advanced Stokes Polarimeter, about a near disc center sunspot region, NOAA AR 9661. We compare the 
characteristics of Stokes profiles in terms of Doppler shifts and asymmetries among the three spectral lines, 
which helps us to better understand the chromospheric lines and the magnetic and flow fields in different 
magnetic regions. The main results are: (1) For penumbral area observed by the photospheric Fe I lines, 
Doppler velocities derived from Stokes / (v\) are very close to those derived from linear polarization profiles 
(yip) but significantly different from those derived from Stokes V profiles (y zc ), which provides direct and strong 
evidence that the penumbral Evershed flows are magnetized and mainly carried by the horizontal magnetic 
component. (2) The rudimentary inverse Evershed effect observed by the Mg I \>2 line provides a qualitative 
evidence on its formation height that is around or just above the temperature minimum region. (3) v zc and v\ p 
in penumbrae and v zc in pores generally approach their i/j observed by the chromospheric Mg I line, which 
is not the case for the photospheric Fe I lines. (4) Outer penumbrae and pores show similar behavior of the 
Stokes V asymmetries that tend to change from positive values in the photosphere (Fe I lines) to negative 
values in the low chromosphere (Mg I line). (5) The Stokes V profiles in plage regions are highly asymmetric 
in the photosphere and more symmetric in the low chromosphere. (6) Strong red shifts and large asymmetries 
are found around the magnetic polarity inversion line within the common penumbra of the 5 spot. We offer 
explanations or speculations to the observed discrepancies between the photospheric and chromospheric lines 
in terms of the three-dimensional structure of the magnetic and velocity fields. This study thus emphasizes the 
importance of spectro-polarimetry using chromospheric lines. 

Subject headings: Sun: activity — Sun: atmospheric motions — Sun: chromosphere — Sun: magnetic fields 
— Sun: photosphere — sunspots — line: profiles — techniques: polarimetric 



1. INTRODUCTION 

Solar photospheric magnetic fields and their associated 
dynamics have been extensively studied. However an 
understanding of their vertical variation from the pho- 
tosphere though the chromosphere to the corona needs 
better comprehension. Direct measurement and diagnosis 
of magnetic and flow fields in higher layers, especially in 
the near force-free chromosphere can play an important 
role in fully understanding their three-dimensional (3D) 
structures. Nevertheless, such efforts were relatively rare due 
to the paucity of chromosp heric spectral lines with s uitable 
Zeeman-split sensitivity (Dalgar no & Layzerl 119871) and 
the complicated dynamics and topology of the magnetized 
plasma in this particular layer. Taking advantage of modern 
observational techniques and recognizing the critical need to 
explore chromospheric magnetic fields, increasing attention 
has recently bee n paid to the spectro-polarimetry of the chro- 
mosp he re (e.g.. | Socas-Navarro et al.l I2000t iZhang & Zhang) 
1 2000b iTruiillo B uenoet alJ l2002t lLopez A riste & C asinil 
2002t ISolanki et alj|2003t iLeka & Metcalj 120031: l Lagg et all 
2004; Balasubramaniam et al. 2004), which reveal the 
height-dependent variation of different types of magnetic 
features. Those analyses were made using magnetically 
sensitive spectral lines formed in the chromospheric levels, 



such as Ca II 849.8 and 854.2 nm, He I 1083.0 nm and D 3 
line, Na I Di 589.6 nm, Ha and H/3 lines. In particular, the 
Mg I b2 517.27 nm line is favorable for spectro-polarimetry 
in the low chromosphere because its core forms in a narrow 
region near or right above the temperature minimum region 
(TMR ) with a relatively large Lande g e /f factor of 1.75 
(e.g., [B riand & Solanki 199 5D. Theoretica l calculation 
( Altro ck & Canfieldl 1 1974 iLitesetalJ [19881 lMauasetal.1 
1988) and observational analy sis (Briand & So lanki 1199 
Briand & Martinez Pilled 1200 U iGosain & Choudharvl [20031 
have been carried out on Mg I b2 Stokes polarization spectra, 
which corroborate that it is well suited for probing the 
magnetic and flow fields in the low chromosphere. 

Doppler shifts of Stokes / and Q,U,V signals can be used 
to investigate mass flows in and aroun d magnetic elements 
dSolankil Il986t ISolanki & Pahlkd [l988). In contrast to the 
Stokes / that represents the total intensity of light, the Stokes 
Q,U (the linearly polarized component) and V (the circularly 
polarized component) have contribution only from the magne- 
tized plasma. Thus Doppler shifts of the Stokes / and Q,U,V 
profiles reflect average line-of-sight (LOS) velocities of plas- 
mas in the whole resolution el ement and tho se within mag- 
netic regions, respectively (e.g.,[Splanki 198g). 

Asymmetries of Stokes profiles are a powerful diagnos- 
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tics of the heig ht and spatial propertie s of magnetic and 
flow fields (e.g., Grigorev & Katz Hl975l) . Under a restric- 
tive Milne-Eddington atmosphere where magnetohydrody- 
namic (MHD) conditions (e.g., plasma velocity, magnetic 
field vector, and temperature) are constant with dept h and 
the source function varies linearly with optical depth (|TJnnol 
ll956tlLandi Deglinnocenti & Landi Deglinnocentil 19851) . the 
Zeeman-split Stokes I,Q,U profiles are symmetric, while V 
is antisymmetric about the line center. Most of the observed 
Stokes spectra, however, deviate from the ideal forms and 
exhibit notable difference between the blue and red wings 
in amplitude and/or area. Therefore the Stokes asymme- 
try is able to reveal the inhomogeneou s nature of the so- 
lar atmosp here (e.g.. Illling et al]|1974a|lb|). Direct obse rva- 
tions (e.g., Balasubraman iam et al.lll997tlSigw arth 20 0lb and 
mode l ing of synthetic Stokes profiles ( e .g.JSolanki & P ahlke 
19881: iSanchez Almeida & Litesl 119921: [S olanki & Montavo'n] 
1993; Sanchez Almeida et al. 1996; Grossmann-Doerth et al. 
2000; Steiner 2000) have found possible physical conditions 
for producing asymmetric Stokes spectra, such as gradients 
of flow velocity and/or magnetic field vector along the LOS 
direction and nonuniform atmosphere that contains two or 
more distinct magnetic and/or flow component s in one reso- 
lution element. A more comprehensive study by Lop ez Aristd 
(2002) based on the analytical solution for the generalized 
transfer equation explicitly demonstrates that the velocity gra- 
dients are the necessary and sufficient condition for Stokes 
profile asymmetries. The addition of other effects or gradi- 
ents of other MHD parameters (e.g., magnetic field and tem- 
perature) could enhance the asymmetries and give different 
characteristics of asymmetries. For example, strong Stokes / 
and V asymmetries are usually found around magnetic flux 
concentrations in the photosphere, which is attributed to the 
prese nce of the "magnetic canopy" (iGrossmann-Doerth et al.1 
Il988ll l989: Solanki 1989; Leka & Steiner|200il). Strong gra- 
dients in both magnetic and flow fields are naturally generated 
because the "magnetic canopy" spatially separates the field- 
free convective downdrafts from the more static flux tube that 
fans out above the photosphere. 

The aforementioned analyses of Stokes profiles predomi- 
nantly use photospheric spectral lines. In retrospect, only few 
such attempts using ch romospheric Stokes spec tra have been 
reported. As examples, Briand & Solanki] (119981) inferred fast 
down flows concentrated in narrow lanes surrounding the flux 
tube below the "magnetic canopy" by using Stokes V asym- 
metries of Mg I b? line in a dditional to photospheric lines. 
iGosain & C houdhary1 (120031) compared the Stokes V asym- 
metries between photospheric Fe I lines and chromospheric 
Mg I lines and found differences in their spatial distribution 
within a simple sunspot. Additional investigation of chromo- 
spheric Stokes spectra in conjunction with photospheric mea- 
surements would improve our understanding of 3D magnetic 
topology and plasma flows. 

In this paper, we study Doppler shifts and asymmetries 
of Stokes profiles of three spectral lines simultaneously ob- 
served in an active region close to the disk center. The three 
spectral lines, Fe I 630.25 nm (g e ff = 2.5), Fe I 630.15 nm 
(geff = 1-67), and Mg I b2, are formed at partially overlap- 
ping atmospheric layers spanning from the photosphere to the 
low chromosphere. Even though both are regarded as photo- 
spheric lines, the Fe I 630.15 nm line often forms above the 
630.25 nm line (e.g., lKhomenko & Colladosll2007l) . Combin- 
ing data from multiple atmospheric heights using these spec- 
tral lines thus enables us to shed new light on the properties of 
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FIG. 1 .— (a) The Fe I 630. 15 nm continuum image of the NOAA AR 9661 
on 2001 October 17. White and black contours outline the positive and neg- 
ative polarities, respectively, (b) The corresponding longitudinal magnetic 
field mimicked by the peak value of Stokes V blue lobe. 

magnetic and flow structures associated with different mag- 
netic features. The plan of this paper is as follows: in §§|2] 
and [3] we introduce the observations and present the proce- 
dure for data reduction. In § |4]the results of comparison of 
profile shifts and asymmetries among the three spectral lines 
are described and discussed. We summarize our major find- 
ings and present some prospects for future studies in §[5j 

2. OBSERVATIONS 

We obtained the spectropolarimetric data sets using the 
76 cm Dunn Solar Telescope and its associated Adaptive Op- 
tics system (iRimmeld 120041) . along with the High Altitude 
Observatory/National Solar Obser vatory Advanced Stokes 
Polarimeter(HAO/NSO ASP. lElmore et^H^9^ lLit^lT996L 
Skumanich et all 1 19971) . The observations were obtained on 
2001 October 17 under excellent seeing conditions. The 
ASP is a slit-grating spectropolarimeter. The 1.6' x 0.6" slit 
was oriented in the solar north-south direction and scanned 
in the east-west direction. Each whole scan took 280 steps 
with a step size of 0.6", which results in a field of view 
(FOV) of 1.6' x 2.8' that covers the near disk center (N14°, 
W6°, [i =cos6> = 0.96) active region NOAA 9661 as shown 
in Figure Q] It contains a S sunspot at left, an a sunspot at 
right, as well as several pores and plages lying in between. 
The obtained spatial sampling is about 0.42" pixel" 1 in the 
north-south direction and 0.6" pixel" 1 in the east-west direc- 
tion. Full Stokes I,Q,U,V spectra were taken simultaneously 
in two spectral bands centered at 630.2 nm and 517.27 nm 
by two ASP cameras. The dispersions for the two spectral 
bands are 1.27 and 1.02 pm pixel" 1 , respectively. Each whole 
scan took about 25 minutes. A total of five whole scans 
were continuously taken within 2.5 hours, during which no 
significant evolution of the magnetic structures was found. 
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Other information about the observatio n run can be found in 
IChou dharv & Balasubra maniaml d2007l) . 

3. DATA REDUCTION 

We applied the standard ASP cal ibration procedures 
dLites et alJl991fclSkumanich et al.ll997l) to the data sets. The 
calibrated Stokes /, Q,U,V spectra were normalized to the 
quiet Sun continuum intensity I c that was obtained by fitting 
the Kitt Peak FTS atlas to the observed Stokes / profiles. The 
data sets from the five whole scans were then integrated to in- 
crease the signal to noise ratio (S/N). Furthermore, we calcu- 
lated the standard deviation a in near continuum wavelength 
range for each integrated Stokes /, Q, U, V profile to represent 
the profile noise level. Only profiles with signal amplitude 
greater than la (an empirical S/N threshold that performs well 
in our data analysis) were used for this study. 

We extracted the parameters for analyzing the shift and 
asymmetry of Stokes profiles using the following procedures. 

( 1 ) Stokes I LOS velocities v, . W e employed a Fourier phase 
method (see Schmid Fet al.ll 1999h to determine the line core 
wavelength. We then converted the wavelength shifts to ve- 
locities (i>i) by using the Doppler formula. As a frame of refer- 
ence, we use the average velocity of a small area in the nearly 
motionless umbra. 

(2) Stokes V LOS velocities v zc . Normalized Stokes V pro- 
files were smoothed to remove local noise. They were then 
divided into two cases: normal profiles that have two opposite 
lobes and one zero-crossing (ZC) point in between, and abnor- 
mal profiles that have no or multiple ZCs within a wavelength 
range corresponding to ±5 kms -1 from Stokes / line center. 
Using this method, there is a small chance that we may im- 
properly classify normal profiles into an abnormal class. For 
exa mple , strong magneto-optical effects (see Fig. [3]: and text 
in § 14. Q produce Stokes V profiles with 3 ZCs, although the 
profile is essentially normal. Such improperly classified cases 
only account for a small portion (< 0.3%) in the data and do 
not affect the analysis. For normal profiles, the ZC shifts are 
converted to LOS velocities (v zc ) in a sa me m anner as Vj. The 
abnormal profiles will be discussed in § 14.11 

(3) Linear polarization (LP) LOS velocities v\ p . The com- 
binati on (Q 2 +U 2 ) x l 2 is a m easure of the overall LP magnitude 
(e.g.. lLeka & Steinerll2001l) . The LP profiles were smoothed 
to remove local noise. It is difficult to measure the position of 
the central 7r component as it is usually small and sometimes 
even vanishes (see Figure [7ji. Instead we measure the peak 
positions of the two a components and use their mid-point 
as the position of the LP profile. We then converted the LP 
profile shifts to LOS velocities (yip) in a same manner as v,. 

(4) Amplitude asymmetry 5a and area asymmetry 8 A. The 
peak amplitudes of the blue (at,) and the red (a r ) lobes of nor- 
mal Stokes V profiles and those of the two a components of 
LP profiles are determined. The areas of blue (Ab) and red 
(A r ) lobes are obtained by integrating over the same wave- 
length range (0.02 nm) on both sides of the ZC wavelength for 
Stokes V profiles or the mid-point of the two a components 
for LP profiles. The integrating range not only covers the ma- 
jority of the Stokes V or LP signal corresponding to Stokes 
/ Doppler core, but also rules out the blends from theneigh- 
boring telluric lines. Following lSolanki & Stenflol (119841) . the 
amplitude asymmetry 5a and area asymmetry 8A are derived 
using the following formulae: 



\A b \+\A r 



4. RESULTS AND DISCUSSIONS 

Figure [2] shows the Stokes /, V and LP images constructed 
by line core intensity and integrated Stokes V and LP sig- 
nals, respectively. Images of the three spectral lines are ar- 
ranged following the order of their formation heights (increas- 
ing height from bottom to top). The orange points are the lo- 
cations of noisy profiles with S/N < 7. The LP spectra have 
more noisy profiles than Stokes V, because the Zeeman sig- 
nals for LP are essentially lower and near the disk center the 
transverse magnetic fields are weaker than the longitudinal 
fields for most of the areas in the low atmosphere. Compared 
to the photospheric Fe I spectra, the chromospheric Mg I po- 
larimetric (Stokes V and LP) spectra generally have lower S/N 
thus more useless profiles. The reasons are mainly due to the 
larger line width and the lower magnetic field strengths in the 
chromosphere. In the following, we first address the spatial 
distribution of abnormal Stokes V profiles, which are impor- 
tant fo r evaluating the complexity of magnetic and flow fields 
(e.g.. lMickevl[T98l ISigwarthll200ll and refer ences therein). 
We then present the results of Doppler shifts and asymmet- 
ric properties of the Stokes profiles, where all the noisy and 
abnormal profiles are excluded from analysis. 

4. 1 . Abnormal Stokes V Profiles 

In the middle column of Figure |2j locations with abnor- 
mal Stokes V profiles are marked using colored points. They 
are divided into three kinds according to the number of ZC 
within a range of ±5 kms" 1 from their Stokes / line cen- 
ters as shown in Figure [3] Profiles with two ZCs (green) are 
most frequently seen, which have three lobes with the cen- 
tral one opposite to the other two (Fig. O?). They concen- 
trate in the regions of magnetic polarity inversion line (MPIL) 
and the outer penumbral boundaries. They are mainly caused 
by the presence of mixed polaritie s associated with differ- 
ent velocities (e.g.. ISigwarthl 2001), whereby the two polar- 
ities can be present beside each other (i.e., within one res- 
olution elejmejit]_OT_aj£nj^ 

regio n dSanchez Almeida & Litesll 992: Sola nki & Montavonl 
1993). In outer penumbral regions, the number of two ZCs 
decreases with height, which implies that the mixed polarity 
effect is weaker at higher atmospheric levels. This is prob- 
ably due to the fact that the orientation of outer penumbral 
fields changes from horizontal or even downward to more ver- 
tical w hen they extend from the photosphere to the chromo- 
sphere (Choudh ary & Balasubraman iam 2007, and references 
therein). However, in the central part of the MPIL of the 8 
spot there are more locations with two ZCs in the Mg I image. 
We speculate that this may point to a stronger mixed polarity 
effect around the highly non-potential MPIL in the chromo- 
sphere. 

The central part of the 5 spot's MPIL is also populated by 
locations without ZC (red). Such locations may involve mag- 
netic components carrying high speed flows, which cause the 
ZC to shift far from the line center (see Fig. (3}?). Since the ZC 
position of Stokes V is susceptible to the prese nce of veloc- 
ity gr adients through the line formation region ( Lope z Aristd 
2002) or distortion of line profiles by noise, it is also possible 
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Stokes I Stokes V LP 




■ noise ■ noise « ZC ■ 2 ZCs « 3 ZCs « noise 

FIG. 2. — Stokes / (line core intensity), Stokes V (integration over Stokes V blue lobe) and LP (integration over LP signal) images in the three spectral lines. 
The sunspot umbrae, penumbrae and pores, as identified on the basis of continuum intensity, are outlined by the yellow contours. The orange points on Stokes /, 
V and LP images indicate the locations where the Stokes /, V and LP profile has a S/N less than 7, respectively. The positions with abnormal Stokes V profiles 
(red, green, and blue for zero, two, and three ZCs, respectively) are also denoted on each Stokes V image accordingly. To enhance the visibility, the size of the 
colored points is larger than the real size of the corresponding image pixel. 



that the profiles at these locations are highly asymmetric or 
distorted. 

There exists another kind of profile with three ZCs (blue). 
They could be caused by residual noise or mixed polarity ef- 
fect. Other than this, the three ZCs profiles within the umbra 
seen in 630.25 nm Stokes V image are most likely due to the 
magneto-o ptical effect (see F ig. |3b) caused by strong Zeeman 
effect (e.g.. lWest & Hagyardlll983l) . 

4.2. Doppler Shifts of Stokes I, V and LP Profiles 

In the upper panels of Figure|U we show Dopplergrams (yi, 
v zc and vip) for each spectral line. They are constructed using 
the LOS velocities derived from the shifts of Stokes /, V and 
LP profiles at each spatial point, respectively. We note the 
followings: (1) For the two photospheric lines, the spatial dis- 
tribution of penumbral Evershed flow in vi Dopplergrams is 
similar to that in vi p but significantly different from v zc Dopp- 
lergrams. (2) The variation of Evershed flow with height can 
be seen by the three spectral lines. Interestingly, the Ever- 
shed flow reverses its sign from the Fe I to Mg I Doppler- 
grams in some places as illustrated by a center-side penumbra 
{white box), which is known as th e inverse Ever shed effect 
observed in the chromosphere (e.g., Solanki 2003). While in 
some other places, the sign reversal has not yet happened or 
is rudimentary. This provides a qualitative evidence on the 
formation height of the Mg I \>2, which is at or just above the 
TMR. (3) The Fe I vi Dopplergrams show that small pores and 
plages lying in between the two sunspots are generally sur- 



rounded by red-shifted downflows, some of which do not ap- 
pear in the Mg I v t Dopplergram as pointed by white arrows on 
the top three panels of Figure|4] In contrast, the Fe I v zc Dopp- 
lergrams show ZC redshifts all over the plage regions. We 
suspect that the red-shifted v zc signal may not represent real 
down flows in flux tubes but could result from the combina- 
tion of the susceptibleness of ZC position to the large Stokes 
V asymmetry often observed by photospheric lines in plage 
regi ons and the limited spectral resolution due t o smoothing 
(see lSolanki & Stenfldfl986tlL6pez Arist3l2002h . 

We quantitatively compare the vi with v zc and v\ p in differ- 
ent magnetic regions for each spectral line by making scatter 
plots of Vi vs. v zc and v-, vs. v\ p shown in the forth and fifth 
row of Figure |U respectively. Since the LP signal is only 
strong enough in penumbrae and MPIL, velocities from these 
regions are plot in Vj vs. v\ p panels. Each group of data points 
are fit with a form of v zc = kv-, + b or v\ p = kvj + b, and their 
linear correlation coefficients (C.C.) are also obtained. We 
discuss the scatter plots for different regions as follows. 

Pores — For the two photospheric Fe I lines, \l> zc \ < |z/;| (i.e., 
k < 1) is observed with moderate correlation (C.C.= 0.56 and 
0.70, respectively). For the Mg I line, \v zc \ « \vj\ (i.e., k rs 1) 
with strong correlation (C.C.= 0.85) is observed. It is well 
established that the magnetic filling factor (MFF) is almost 1 
and the magnetic fields are nearly vertical in the pores. This 
is even true for pores in heights above the photosphere, which 
explains that v zc show large similarity with Vj as seen by the 
Mg I line. While down to the photosphere near the edge of 
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FIG. 3. — The sample abnormal Stokes V profiles defined in this study with zero, two, and three ZCs within ±5 kms from their Stokes / line centers indicated 
by the vertical solid line. The vertical dotted lines mark the wavelength range corresponding to ±5 kms~' from Stokes / line center. 



pores the field lines may be more inclined and the MFF may 
be less than 1 in such an observatio n with median spatial res- 
olution (e.g., Leka & Steiner 2001). Strong and narrow field- 
free down flows with speed larger than that within the pores 
have been found immediately around por es in the photosphere 
(Sankarasubramanian & Rimmele 2003), which gives rise to 
| Uj | instead of | v zc \ . This might be responsible for the observed 
\v zc \ < | v; | in the Fe I lines. 

Penumbrae — The penumbral magnetic fields are known to 
consist of a more vertical component and a more horizon- 
tal compo nent in the photosph ere (interlocking-comb struc- 
ture, e.g., Ilchimoto et al .1120071) . For this region close to the 
disk center, v zc have more contribution from the more verti- 
cal component than from the other component. On the con- 
trary, vi p have more contribution from the more horizontal 
component. Figure [4] shows that v, has much better spatial 
and magnitude correlation with v\ p than with v zc for the two 
Fe I lines. Moreover, the slopes of Vj vs. v\ p (k= 0.87 and 
0.86, respectively) are close to 1 and much larger than those 
of i>i vs. v zc (k= 0.33 and 0.35, respectively). In other words, 
the vi p is very close to v-, and their magnitude is much larger 
than \v zc \. These results directly support that the penumbral 
Evershed flows are not only magnetized but mainly carried 
by the horizontal magnetic component as previously found by 
othe r authors using different observations and meth ods (see, 
e.g.. iRimmele & Marinoll2006l: Ilchimoto et al.ll2008[ and ref- 
erences therein). While for the Mg I line, v-, is observed sim- 
ilar to both of Vi p and v zc (i.e., k ps 1) with small correla- 
tion coefficients (C.C.= 0.21 and 0.40, respectively). This 
implies that rising up to the TMR the penumbral magnetic 
and flow structures are different from those in the photo- 
sphere. The large scatter of the data points resulted from nois- 
ier and weaker polarimetric signal of the Mg I line hinders us 
from drawing solid conclusion. We thus speculate that the 
interlocking-comb structure of the penumbral magnetic and 
flow field might be a shallow phenomenon mainly in the pho- 
tosphere. 

MPIL — For a compact elongated region around the MPIL 
and within the common penumbra of the 5 spot, apart from 
noise dominating the Fe I 630.25 v-, and v\ p signal, the Dopp- 
lergrams show strong redshift with a mean speed conspicu- 



ously larger than that of simultaneously observed Evershed 
flows. v zc and v\ p are generally well correlated with v t . The 
linear fittings show that \v\ p \ and \v-\ are close to each other 
and much smaller than \v zc \. The large \v zc \ may be ar- 
tificial signal resulted from large Stokes V asymmetry and 
limited spectral resolution as discussed above. It is also 
possible that the mixed polarities there contain a compo- 
ne nt possessing a high speed downflow, as observed earlier 
bv lMartinez Pillet et al.l (U994). Detailed examination of the 
large redshift requires two-component Stokes inversion of in- 
dividual profiles, which is out of the scope of this study. 

4.3. Asymmetry of Stokes Profiles 

Asymmetries of both Stokes V and LP profiles were cal- 
culated. We present results from Stokes V only because the 
noisier and weaker LP signals produce large uncertainty in de- 
termined asymmetries although their spatial distributions are 
qualitatively similar to those of Stokes V. Figure[5]shows the 
spatial distribution and spectral variation of the Stokes V am- 
plitude asymmetry {8a) and area asymmetry (8A). We also 
plot their histograms for different magnetic regions in Fig- 
ure [6] in order to quantitatively compare the Stokes V asym- 
metries among the three spectral lines. The typical Stokes I,V 
and LP profiles in those regions are also plotted in Figure [7] 
We again discuss various magnetic features separately as fol- 
lows. 

For plages, the red lobes are broader and shallower than 
the blue lobes for the two photospheric Fe I lines, which 
makes highly asymmetric Stokes V profiles thus substan- 
tial positive 6 a and 5 A with 8 a > 5A generally, as pre- 
viousl y found (IStenflo et al.lll984t ISolanki & Stenflol 11984 
[1981 iGrossmann-Doerth et all 1 19961: ISigwarthll200lT) . It is 
believed that the "magnetic canopy" effect and the small 
MFF in photospheric plages play importa nt roles in pro- 
ducing such large positive asymmetries (e.g.. lLeka & Steiner] 
120011) . A noteworthy finding here is that the chromospheric 
Mg I Stokes V profiles are much more symmetric as can be 
seen that th e most probable 8a and 8 A are very small posi- 
tive values. iBriand & Solankil (119981) also observed smaller 
Stokes V asymmetry in the Mg I b2 line compared to a photo- 
spheric line in network regions. It was found that for small or 
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FIG. 4. — Upper panels: Dopplergrams made using Stokes / line center shifts (V, ), Stokes V ZC shifts (y zc ) and LP profiles shifts (vi p ) for the three spectral 
lines. The numbers in unit of kms~'are the mean velocities within the strong and compact redshift region around the MPIL of the <5 spot, and those in the 
penumbral region outlined by the white box. The spatial points with noisy and abnormal profiles are assigned a value of 0. Note that the methods of determining 
Vj, u zc and ui p may break down in some umbral areas for fully split lines. Lower panels: Scatter plots of Vj vs. v- c (forth row) and Vj vs. V] p (fifth row) for 
each spectral line. The color-coded data points stand for different magnetic regions. The slope k and the linear correlation coefficient (C.C.) from the line fits 
u zc = kvi + h or vi p — ki/j + b for each group of data are labeled. 
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FIG. 5. — Images constructed with values of Stokes V amplitude asymmetries (Sa, first row) and area asymmetries (<5A, second row) at each spatial point for 
the three spectral lines. The magnitude of asymmetry is indicated by the scales at right. The spatial points with noisy and abnormal profiles are assigned a value 
ofO. 




FIG. 6. — The distribution of Stokes V amplitude asymmetries (5a, first row) and area asymmetries (<5A, second row) in different magnetic regions for the three 
spectral lines. The black, gray, and dotted lines are for the Fe I 630.25, Fe 1 630.15, and Mg 1 517.27 observations, respectively. 



weak magnetic structures, such as plage and network, the de- 
gree of Stokes asymmetries decreases if MFF increases (e.g., 
Zaver et alj[l990t iMartmez Pillet et all 1 19971: ISigwarth et al.1 
1999). The dependence of amplitude of Stokes V asymme- 
tries on some line parameters, such as Doppler shift, Zee- 
man shift, line strength and width, has also been heuristically 
conside red within the frame of a simpl ified magnetic canopy 
model (Grossmann-Doer th et al.|[l9 89). There the asymme- 
try is expected to be small for a same downdraft velocity if 
the line is broad and not strongly Zeeman split. Our result 
confirms these trends from observational point of view. To 
sum up, both the broader line width and the larger MFF in the 
formation layer of the Mg I b2 line are responsible for the ob- 
served small asymmetry. Together with the aforementioned 
finding that some of the conspicuous downflows surrounding 
the plages disappear in the Mg I v, Dopplergram, we there- 
fore suspect that the "magnetic canopy" of plages could be a 
shallow phenomenon mainly below the TMR, which however 
needs more solid observational evidence to justify. 



For penumbrae (especially their outer parts) and pores, 5a 
and SA are primarily positive and negative in the photosphere 
and in the low chromosphere, respectively. Previous studies 
have shown that penumbrae and pores share some common 
properties. They both present an interface with the quiet sun 
and are surr ounded by vortex convective flows with similar 
patterns (e.g. lDeng et alj |2007). This study hence adds to the 
evidence of another common property, i.e., Stokes V asymme- 
tries and their variation with height. Unlike photospheric lines 
that can be easily synthesized under LTE conditions, the phys- 
ical origins of the asymmetries of chromospheric lines formed 
in Non-LTE atmosphere are far from fully studied. Here we 
can only present a speculation for the sign reversal of Stokes V 
asymmetries based on a logical analog of chromospheric line 
asymmetries to the better known photospheric line asymme- 
tries that are mainly caused by velocity gradients. The change 
in sign of Evershed effect for penumbrae may be responsible 
for some but not all areas, because the Stokes V asymme- 
tries reverse sign all over the outer penumbral borders while 



8 



DENG ET AL. 



Fel 630.15 & 630.25 nm line pair 
plage CS penum. LS penum. pore umbra MPIL 




630.15 630.20 630.25 630.15 630.20 630.25 630.15 630.20 630.25 630.15 630.20 630.25 630.15 630.20 630.25 630.15 630.20 630.25 



Wavelength [nm] 



Mgl 517.27 nm 

plage CS penum. LS penum. pore umbra MPIL 




517.20 517.25 517.30 517.35517.20 517.25 517.30 517.35517.20 517.25 517.30 517.35517.20 517.25 517.30 517.35517.20 517.25 517.30 517.35517.20 517.25 517.30 517.35 



Wavelength [nm] 

FIG. 7. — Typical Stokes /, V and LP profiles of the Fe I line pair and the Mg I b2 line averaged over different magnetic regions: plage, center-side (CS) 
penumbra, limb-side (LS) penumbra, pore, umbra, and the compact strong red-shifted region around the MPIL of the <5 spot. 
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the inverse of Evershed flow only occurs in some penumbral 
sectors. We speculate that the sign reversal of Stokes V asym- 
metries in both outer penumbrae and pores could indicate that 
the velocity in and around them may change drastically from 
the photosphere to the low chromosphere. Figure H] provides 
clues to some of the changes because some downflows in and 
around pores do not hold in the Mg I v-, Dopplergram instead 
more upflows ca n be seen (e.g., in the scatter plot of the Mg I 
line). Moreover, Sankarasubramani an & Rim mele (2003]) ob- 
served pores in different wavelengths and found that veloci- 
ties in the narrow downflow region around pores tend to de- 
crease with height at first, and then change to upflow in the 
low chromosphere. Such tendenc y is also visible in the MHD 
simulation of pores (see Fig. 3 of iLeka & Ste iner 200 1]). The 
Stokes V asymmetries change from positive for the Fe I lines 
to negative for the Mg I line, if indeed all line asymmetries 
have a velocity gradient origin, thus provides important clues 
on the 3D structure of the flow field in and around pores and 
outer penumbrae. 

For umbrae, since they lack of mass flows and have MFF of 
~1 in the low photosphere, both of the most probable Sa and 
SA of the Fe 1 630.25 nm line are close to 0. The large Zeeman 
splitting of the line in stron g field may also contribute to th e 
observed small asymmetry (Grossmann-Doerth et al. 1989). 
However, the distributions of Sa and SA of the other two spec- 
tral lines show a clear shift toward negative. Besides the noise 
induced inaccuracy of the determined asymmetry due to low 
light level in the umbra, it is also possible that as height in- 
creases umbrae become less uniform and static compared with 
their structure in the low photosphere. 

For the compact region with strong redshifts around the 
MPIL of the S spot, both da and SA have prominent values that 
increase with height. It seems from Figure[7]that the red wings 
of the Stokes /, V, and LP profiles of the two photospheric 
lines obviously have further redward extension. Meanwhile, 
the blue damping wings of the Mg I 517.27 nm line are sur- 
prisingly enhanced, which remains a puzzling issue. We con- 
jecture that the observed large asymmetries could be due to 
the combination of the following two effects: the mixed po- 
larities or flows, and the large gradients in velocity, magnetic 
field vector, or temperature near the MPIL region. Consider- 
ing all the unusual properties observed in this region, we sug- 
gest that the magnetic and flow fields around the MPIL of the 
S spot are very complex, which deserve a deeper investigation 
in the future. 

5. SUMMARY 

We have presented a comprehensive study of the Doppler 
shifts and asymmetric properties of the Stokes profiles, the 
uniqueness of which is especially in the combination of mag- 
netic sensitive spectral lines that form at different heights 
spanning from the photosphere to the low chromosphere. 
Our quantitative analysis benefits from the simultaneous and 
high resolution spectro-polarimetric observation that include 
a chromospheric line. We summarize the major findings and 
interpretations as follows. 

1 . Vj is very close to v\ p but significantly different from v lc 
in penumbral areas near disc center for the photospheric 
Fe I lines, which provides direct and strong evidence 
that the penumbral Evershed flows are magnetized and 
mainly carried by the horizontal magnetic component. 

2. The rudimentary inverse Evershed effect observed by 



the Mg I b2 line provides a qualitative evidence on its 
formation height that is around or just above the TMR. 

3. v zc and vi p in penumbrae and v zc in pores generally ap- 
proach their Vi observed by the chromospheric Mg I 
line, which is not the case for the photospheric Fe I 
lines. We speculate that the interlocking-comb struc- 
ture of the penumbral magnetic and flow field might be 
a shallow phenomenon mainly in the photosphere. 

4. The Stokes V asymmetries exhibit similar behavior in 
the regions of outer penumbrae and pores. Both 5a and 
8A evolve from primarily positive values observed by 
photospheric Fe I lines to primarily negative values ob- 
served by the chromospheric Mg I line. 

5. In plage regions where the "magnetic canopy" effect 
is present, Stokes V profiles are highly asymmetric in 
the photosphere but more symmetric in the low chro- 
mosphere. Moreover, some of the photospheric down- 
flows surrounding the plages disappear in the Mg I 
Dopplergram. These suggest that the plage "magnetic 
canopy" could be a shallow phenomenon mainly below 
the TMR. 

6. All the Vi, v zc and vi p Dopplergrams reveal strong red- 
shifts in a compact region around the MPIL and within 
the common penumbra of the S spot. Large posi- 
tive Stokes V asymmetries are also present and even 
strengthen with height. These unusual phenomena are 
manifestations of the complex nature of the 3D topol- 
ogy of the highly non-potential region. 

We conclude that analysis of Stokes profiles using both 
the photospheric and chromospheric spectral lines can fur- 
nish crucial information on the 3D structure of the magnetic 
and flow fields. Similar study using other chromospheric lines 
will be essential to see if the behavior and trend of velocities 
and asymmetries observed by the Mg I bo line is general for all 
the chromospheric lines. The speculations made in this study, 
such as the velocity gradient origin of the chromospheric line 
asymmetries and the thinness of the interlocking-comb struc- 
ture of penumbral magnetic and flow fields, need further in- 
vestigation to justify. To advance firmly towards a sound and 
reliable measurement of magnetic fields in the chromosphere, 
not only more polarimetric observations using chromospheric 
lines are needed, but also a better understanding of their com- 
plicated profile formation conditions is critical. 

We are truly grateful to Dr. S. Solanki for careful reading 
of the paper and many constructive suggestions throughout 
this study. Thanks are also given to Dr. C. Liu and A. Cook- 
son for polishing and improving the manuscript and Drs. L. 
Bellot Rubio, K. E., Rangarajan, and H. Wang for helpful dis- 
cussions. We thank the referee for valuable comments that 
help us to improve the paper. N.D. and D.P.C. were sup- 
ported by NASA grant NNX08AQ32G and NSF grant ATM 
05-48260. We acknowledge the use of the HAO/NSO Ad- 
vanced Stokes Polarimeter. The National Solar Observatory 
is operated by the Association of Universities for Research in 
Astronomy under a cooperative agreement with the National 
Science Foundation, for the benefit of the astronomical com- 
munity. This research has made use of NASA's Astrophysics 
Data System Bibliographic Services. 



10 



DENG ET AL. 



REFERENCES 



Altrock, R. C, & Canfield, R. C. 1974, ApJ, 194, 733 

Balasubramaniam, K. S., Christopoulou, E. B., & Uitenbroek, H. 2004, ApJ, 
606, 1233 

Balasubramaniam, K. S., Keil, S. L., & Tomczyk, S. 1997, ApJ, 482, 1065 

Briand, C, & Martinez Pillet, V. 2001, in Astronomical Society of the 
Pacific Conference Series, Vol. 236, Advanced Solar Polarimetry - 
Theory, Observation, and Instrumentation, ed. M. Sigwarth, 565-+ 

Briand, C, & Solanki, S. K. 1995, A&A, 299, 596 

— . 1998, A&A, 330, 1160 

Choudhary, D. P., & Balasubramaniam, K. S. 2007, ApJ, 664, 1228 
Dalgarno, A., & Layzer, D., eds. 1987, Spectroscopy of astrophysical 
plasmas 

Deng, N., Choudhary, D. P., Tritschler, A., Denker, C, Liu, C, & Wang, H. 

2007, ApJ, 671, 1013 

Elmore, D. E, Lites, B. W., Tomczyk, S., Skumanich, A. P., Dunn, R. B., 
Schuenke, J. A., Streander, K. V., Leach, T. W., Chambellan, C. W., & 
Hull, H. K. 1992, in Presented at the Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference, Vol. 1746, Polarization 
analysis and measurement; Proceedings of the Meeting, San Diego, CA, 
July 19-21, 1992 (A93-33401 12-35), p. 22-33., ed. D. H. Goldstein & 
R. A. Chipman, 22-33 

Evershed, J. 1909, MNRAS, 69, 454 

Gosain, S., & Choudhary, D. 2003, Sol. Phys., 217, 119 

Grigorev, V. M., & Katz, J. M. 1975, Sol. Phys., 42, 21 

Grossmann-Doerth, U., Keller, C. U., & Schuessler, M. 1996, A&A, 315, 
610 

Grossmann-Doerth, U., Schuessler, M., & Solanki, S. K. 1988, A&A, 206, 
L37 

— . 1989, A&A, 221,338 

Grossmann-Doerth, U., Schiissler, M., Sigwarth, M., & Steiner, O. 2000, 
A&A, 357, 351 

Ichimoto, K., Shine, R. A., Lites, B., Kubo, M., Shimizu, T, Suematsu, Y., 
Tsuneta, S., Katsukawa, Y, Tarbell, T. D., Title, A. M., Nagata, S., 
Yokoyama, T., & Shimojo, M. 2007, PASJ, 59, 593 

Ichimoto, K., Tsuneta, S., Suematsu, Y, Katsukawa, Y, Shimizu, T., Lites, 
B. W., Kubo, M., Tarbell, T. D., Shine, R. A., Title, A. M., & Nagata, S. 

2008, A&A, 481,L9 

Illing, R. M. E., Landman, D. A., & Mickey, D. L. 1974a, A&A, 35, 327 
— . 1974b, A&A, 37, 97 

Khomenko, E., & Collados, M. 2007, ApJ, 659, 1726 
Lagg, A., Woch, J., Krupp, N, & Solanki, S. K. 2004, A&A, 414, 1 109 
Landi Deglinnocenti, E., & Landi Deglinnocenti, M. 1985, Sol. Phys., 97, 
239 

Leka, K. D., & Metcalf, T. R. 2003, Sol. Phys., 212, 361 
Leka, K. D., & Steiner, O. 2001, ApJ, 552, 354 
Lites, B. W. 1996, Sol. Phys., 163, 223 



Lites, B. W., Elmore, D., Murphy, G., Skumanich, A., Tomczyk, S., & 
Dunn, R. B. 1991, in 11. National Solar Observatory / Sacramento Peak 
Summer Workshop: Solar polarimetry, p. 3 - 15, 3-15 

Lites, B. W., Skumanich, A., Rees, D. E., & Murphy, G. A. 1988, ApJ, 330, 
493 

Lopez Ariste, A. 2002, ApJ, 564, 379 
Lopez Ariste, A., & Casini, R. 2002, ApJ, 575, 529 
Martinez Pillet, V., Lites, B. W., & Skumanich, A. 1997, ApJ, 474, 810 
Martinez Pillet, V., Lites, B. W., Skumanich, A., & Degenhardt, D. 1994, 
ApJ, 425, LI 13 

Mauas, P. J., Avrett, E. H., & Loeser, R. 1988, ApJ, 330, 1008 

Mickey, D. L. 1985, Sol. Phys., 97, 223 

Rimmele, T, & Marino, J. 2006, ApJ, 646, 593 

Rimmele, T. R. 2004, in Presented at the Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference, Vol. 5490, Advancements 
in Adaptive Optics. Edited by Domenico B. Calia, Brent L. Ellerbroek, 
and Roberto Ragazzoni. Proceedings of the SPIE, Volume 5490, pp. 34-46 
(2004)., ed. D. Bonaccini Calia, B. L. Ellerbroek, & R. Ragazzoni, 34^16 

Sanchez Almeida, J., Landi degl'Innocenti, E., Martinez Pillet, V, & Lites, 
B. W. 1996, ApJ, 466, 537 

Sanchez Almeida, J., & Lites, B. W. 1992, ApJ, 398, 359 

Sankarasubramanian, K, & Rimmele, T. 2003, ApJ, 598, 689 

Schmidt, W., Stix, M., & Wohl, H. 1999, A&A, 346, 633 

Sigwarth, M. 2001, ApJ, 563, 1031 

Sigwarth, M., Balasubramaniam, K. S., Knolker, M., & Schmidt, W. 1999, 
A&A, 349, 941 

Skumanich, A., Lites, B. W., Martinez Pillet, V, & Seagraves, P. 1997, 
ApJS, 110, 357 

Socas-Navarro, H., Trujillo Bueno, J., & Ruiz Cobo, B. 2000, Science, 288, 
1396 

Solanki, S. K. 1986, A&A, 168, 311 

— . 1989, A&A, 224, 225 

— . 2003, Astron. Astrophys. Rev. 11, 153 

Solanki, S. K, Lagg, A., Woch, J., Krupp, N., & Collados, M. 2003, Nature, 
425, 692 

Solanki, S. K., & Montavon, C. A. P. 1993, A&A, 275, 283 
Solanki, S. K, & Pahlke, K. D. 1988, A&A, 201, 143 
Solanki, S. K, & Stenfio, J. O. 1984, A&A, 140, 185 
— . 1985, A&A, 148, 123 

— . 1986, A&A, 170, 120 

Steiner, O. 2000, Sol. Phys., 196, 245 

Stenfio, J. O., Solanki, S., Harvey, J. W., & Brault, J. W. 1984, A&A, 131, 
333 

Trujillo Bueno, J., Landi Degl'Innocenti, E., Collados, M., Merenda, L., & 

Manso Sainz, R. 2002, Nature, 415, 403 
Unno, W. 1956, PASJ, 8, 108 

West, E. A., & Hagyard, M. J. 1983, Sol. Phys., 88, 51 

Zayer, I., Stenfio, J. O., Keller, C. U., & Solanki, S. K. 1990, A&A, 239, 356 

Zhang, M., & Zhang, H. 2000, Sol. Phys., 194, 19 



